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Abstract 

Using covariance analysis, we quantify the correlations between the interac¬ 
tion parameters in a transport model and the observables commonly used to 
extract information of the Equation of State of Asymmetric Nuclear Matter 
in experiments. By simulating 124 Sn+ 124 Sn, 124 Sn+ 112 Sn and 112 Sn+ 112 Sn re¬ 
actions at beam energies of 50 and 120 MeV per nucleon, we have identified 
that the nucleon effective mass splitting are most strongly correlated to the 
neutrons and protons yield ratios with high kinetic energy from central col¬ 
lisions especially at high incident energy. The best observable to determine 
the slope of the symmetry energy, L, at saturation density is the isospin diffu¬ 
sion observable even though the correlation is not very strong (~0.7). Similar 
magnitude of correlation but opposite in sign exists for isospin diffusion and 
nucleon isoscalar effective mass. At 120 MeV/u, the effective mass splitting 
and the isoscalar effective mass also have opposite correlation for the dou¬ 
ble n/p and isoscaling p/p yield ratios. By combining data and simulations 
at different beam energies, it should be possible to place constraints on the 
slope of symmetry energy (L) and effective mass splitting with reasonable 
uncertainties. 

Keywords: nucleon effective mass and effective mass splitting, symmetry 
energy, heavy ion collisions,covariance analysis 


Knowledge about the isospin asymmetric nuclear equation of state (EoS) 
is of fundamental importance to our understanding of nature’s most asym- 
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metric objects including neutron stars and heavy nuclei composed of very 
different numbers of neutrons and protons. Theoretically, there are two 
microscopic approaches to describe the EoS of nuclear matter. One ap¬ 
proach starts from a realistic two-body free NN interactions [l, 0] as input to 
the relativistic Dirac-Bruckner-Hartree-Fock (DBHF) and its nonrelativis- 
tic counterpart Bruckner-Hartree-Fock (BHF)|j], 3 , 1 , 0 ] and chiral ef¬ 

fective field theoryj8j. Another one is to use effective density-dependent 
many-body interactions such as the zero-range Skyrme interaction|9l llOL 
finite-range Gogny interaction jl2j] and effective Lagrangian [ll3[ flil . TlaL 
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as inputs leading to Skyrme-Hartree-Fock (SHF) 


Gogny- H art ree- 


Fock 17] and Relativistic-Hartree-Fock (RHF) approaches [131. 14.1151 Il6 


all interactions, the effective Skyrme interaction are more commonly used in 
nuclear structure, reactions and astrophysics studies as the effective Skyrme 
interactions are relatively simple mathematically to make it computationally 
feasible fl8|] and contain sufficient physics to allow quantitative description 
of heavy nuclei. Furthermore, the Fock term in the non-relativistic SHF 
approach can be managed computationally while the Fock terms in Rela¬ 
tivistic HF approaches are difficult to compute. The interactions parameters 
are usually obtained by fitting the properties of symmetric nuclear matter 
(such as the saturation density and the corresponding energy per nucleon 
and its incompressibility), properties of asymmetric nuclear matter (such as 
symmetry energy and isovector effective mass at saturation density), finite 
nuclei properties (such as binding energies and r.m.s radius of selected set of 
doubly magic nuclei) etc. Q. In Ref jlij ]. 240 Skyrme parameter sets that 
fit the ground-state properties of stable nuclei, symmetric and asymmetric 
nuclear matter are compiled. The large number of parameterization arises 
in part because there are strong correlations between individual parameters 
or groups of parameters, that fit particular physical properties of the many- 
body nuclear system. These sets lead to very different Equation of States of 


pure neutron matter 20, 2l|, 22] which may have different incompressibility 


K 0 = 9symmetry energy coefficient S 0 = S(p 0 ), slope of symmetry 


3 Po 9S q ^ I pen isoscalar effective mass m , 


energy L = 
and isovector effective mass m* = 


m _ / -| | 2m a E 

?77* k h 2 dr A) \P 0 
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-2—, where k is the enhancement factor 

1+Rr 2 - 


24]. In this work, the magnitude of 


of the Thomas-Reich-Kuhn sum rule 
effective mass splitting, (m* — m*)/m was not used as input variables since 
its form is much more complicated to be incorporated into the code. In¬ 
stead we use fj = ±( Jr - jjjjr) = ^ where 5 = (p n - p p )/(p n + p P ), 
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n , m* and m are the neu- 


p n and p p are the neutron and proton density, m ; 
tron, proton effective mass and free nucleon mass. In Skyrme Hartree-Fock 
approximation [To], [□], |l9|, fj increases with increasing density, but is in¬ 
dependent of the momentum and kinetic energy. In the DBHF and RHF 

fi not only depends on the density but also on 


approximations [161 . 
the kinetic energy of the in-medium nucleon. 

Some properties of symmetric nuclear matter, such as K 0 = 230±30MeV 17|, 
19, 27] and m*/m = 0.65 — 0.9 jT3, 19, 23, 28, 0, 3o|, have been extracted 


from isoscalar collective vibrations, giant quadrupole resonance and heavy 
ion collisions measurements. In addition, constraints on E 0 (p ) and pressure 
P have been obtained in density regions ranging from saturation density to 
five times normal densities using collective flow and kaon production data 
in energetic nucleus-nucleus collisions 31, [32, [33]. To obtain information of 


the symmetry energy with heavy ion collision data, the symmetry potential 
used in transport models is changed by varying its input parameters, corre¬ 
sponding to the different values of So and/or L in the expression of density 
dependence of symmetry energy. The results of the calculations are then 
compared with data to find the best parameter sets. Recently, a consistent 
set of constraints on the symmetry energy near saturation density between 
S 0 , and its slope, L, has been obtained from observables measured in both 
nuclear structure and nuclear reaction experiments (iil. lilil . l35i 36 . 


In Skyrme parametrizations, the symmetry energy are correlated to other 
parameters such as those associated with the nucleon effective mass m* and 
isovector effective mass m*. Such correlations would affect the uncertainties 
of symmetry energy constraints obtained from heavy ion collision data. In 
order to achieve the goal of obtaining precise and accurate symmetry energy 
constraints, one has to identify what experimental observables are crucial 
for better constraining the interested physical quantities in the theoretical 
models [37]. Ideally, one would do a global chi-square analysis using existing 
data to obtain the best set of model parameters. Then a covariance matrix 


can be obtained between any two model parameters using a chi-square fit 38 


Currently, this is not feasible considering the intensive CPU time needed to do 
transport model calculations. As a start to tackle this issue, we propose to use 
12 parameters sets to perform covariance analysis to quantitatively examine 
the correlations between model parameters A and observables B commonly 
used in experiments. The linear-correlation coefficient Cab between variable 
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( 1 ) 


A and observable B is 
Cab 
cov(A , B) 


a(X) 


< X > 


calculated as follows (38J: 

cov(A 1 B ) 

<j{A)<j{B) 

i 

= jrT, x *’ i = 1 ’ N - 


( 2 ) 

( 3 ) 

( 4 ) 


cov(A,B) is the covariance, <j(X) is the variance. Cab = ±1 means there is 
a linear dependence between A and B } and Cab = 0 means no correlations. 

We use the Improved Quantum Molecular Dynamic Model which incor¬ 
porates the effective Skyrme interactions (ImQMD-Sky) 39] to simulate the 
collisions of heavy ions with parameter sets listed in Table I. Ai represents 
the ith parameter set of the transport variable A where A=A'o, So, L, m*, 
or fj used as input to the ImQMD-Sky. In ImQMD-Sky, the nucleonic 


m 


potential energy density is Ui oc + u m d , where 
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a p~ 
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P P 


,r/+1 
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2 p 0 V + 1 Po 2p 0 

Qsur,iso f ^-7 f _ _ ^-|2 


(Vp) 2 


Po 
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+A sym p 2 6 2 + B sym pv +1 5 2 


( 5 ) 


and the energy density of Skyrme-type momentum dependent interaction are 
written based on its interaction form 5(r\ — r 2 )(pi — P 2 ) 2 0,0, 


U"md ^md^Pn^n) ^md^Pp^~p) ( 6 ) 

= C 0 J (PpcPpf (f , p)f{r , p) (p- p) 2 + 

Do j d :y pd :i p' [f n (f, p) f n (f, f) (p -ff 
+f P if, p) f p (r. , f) (p - fj 2 ]. 

The 9 parameters a, /3, 77 , A sym , B sym , C 0 , D 0 , g sur , g sur ,iso used in ImQMD- 
Sky can be derived from standard Skyrme parameter sets with 9 parameters 
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{to, h, t 2 , t 3 , X 0 , Xi, x 2 , X 3 , cr} (39!, Bq 
terms are set as g sur = 24.5 MeVfm 2 and g 


The coefficients of the surface 

_ so = —4.99 MeVfm 2 which 

are the same values derived from SLy4 parameter set 0- Varying g sur and 
gsur,iso for different Skyrme interactions have negligible effects on the ex¬ 
perimental observables at intermediate energy. By using the relationship 
which was derived in reference (4ll. 0], the reduced 7 Skyrme parameter sets 
{a, (3, T), A sym , B sym , C 0 , Dq} can be replaced by the parameter sets {p 0 , 
Eq, K 0 , Sq, L, m*, m*} which is directly related to the properties of nu¬ 
clear matter at saturation density po- Choosing the experimental values of 
Po = 0.16/m -3 , E 0 = —16 MeV, the parameter sets are further reduced to 
5 variables, A=A' 0 ,S'o,L,m*,m*. As explained above, to simplify the coding, 
the input variables used in ImQMD-Sky are A —K 0 ,S 0 ,L,m*,fi(m*,ml). 

For HIC observables, we adopt the ratios constructed from nucleon spec¬ 
tra. Most transport models cannot describe accurately the absolute yield 
of free nucleons due to models inadequacies in describing light clusters ( 43 I. 
44,145]. This problem can be largely alleviated by constructing ’’coalescence 
invariant” quantities, i.e., nucleon observables summed over all light clusters, 


which show much better agreement between theory and experiment [36l. 146 


471 [48]. We construct the coalescence invariant (Cl) nucleon yield spectra 


and their ratios in the same way as in Ref|ii I 47 L [iq by combining the nu 
cleons in the light particles and free nucleons at given kinetic energy per 
nucleon as follows, 


dYj(n ) 
dE cm JA 
dYj(p) 
dE c m jA 


Z. N 

E z 


dMi(N, Z ) 
dE c m / A 
dMi(N, Z ) 
dE c rn / A 


(7) 

( 8 ) 


The summation is up to Z < 6 and A < 16 particles in the calculations, 
mu }tjphdty of fragments with neutron number N and proton 
number Z at certain kinetic energy. Such Cl neutron and proton yields 
obtained from ImQMD-Skycalculations with selected Skyrme parameter sets 
reproduce the Sn+Sn data[48[ reasonably well especially at high beam energy. 
For simplicity, integrated Cl neutron and proton yields from reaction i are 
represented as V(n) and Y % (p) respectively in the following. By convention, 
the more neutron-rich reaction is represented by the subscript ”2” in this 
work, i=l,2 represent the reactions 112 Sn+ 112 Sn and 124 Sn+ 124 Sn. 

We simulate 10,000 events for each of the three systems, 124 Sn+ 124 Sn, 
124 Sn+ 112 Sn and 112 Sn+ 112 Sn. Four yield ratios are constructed by inte- 
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grating the Cl nucleon energy spectra with E cm /A > 40 MeV. High en¬ 
ergy nucleons are less influenced by sequential decay and are more sen¬ 
sitive to the effective mass splitting[[39, 49, 50, [5l|, 52, 53]. The angular 
cut, 70° < 6 c .m. < 110°, is also imposed to reduce the contributions to 
the Cl nucleon spectra from reaction residues. The four ratio observables 
are: 1.) single n/p ratio QVR 2 {n /p)=Y 2 {;n)/Y 2 (j)), 2.) double n/p ratios CI- 
DR(n/p)=CI-R 2 (n/p)/CI-Ri(n/p)=CI-R 2 i(n/n)/CI-R 2 i(p/p), 3.) isoscaling 
ratios CI-R 2 i(n/n)=Y 2 (n)/Yi(n) and 4.) CI-R 2 i(p/p)=Y 2 (p)/Yi(p) ratios. 
The experimental isoscaling ratios and the double ratios have the advantage 
that they are minimally affected by detector systematic uncertainties and se¬ 


quential decays[54], [55]. 5.) the fifth and last isospin observable we examine 
is the isospin transport ratiosj56j, Rdiff , that quantify diffusion of the nucle¬ 
ons in the neck region during the nuclear collisions. Rdiff have been used to 
constrain the symmetry energy at subsaturation density. These constraints 
are consistent with constraints obtained in nuclear structure experiments ( 22 | . 

The isospin diffusion observable, Rdiff is defined as [56j, Rdiff = (2Y — 
X aa — X bb )/(X aa — X bb ) where X = X ab is an isospin observable. In this 
work, we use the subscripts a and b to denote the projectile (first index) and 
target (second index) combination. By convention, aa and bb represent the 
neutron-rich and neutron-poor reactions, respectively. In theoretical calcula¬ 
tions, X is the isospin asymmetry of the emitting source, which is constructed 
from all the emitted nucleons and fragments with velocity cut v > 0-5v b ™ m . 
It is assumed to be linearly related to the experimental observable X exp as 
discussed in [56|, 57]. Rdiff are constructed with at least three reaction sys¬ 
tems to cancel the drift and retain the information of the diffusion. Since 
the changes of isospin asymmetry of emitted source mainly come from the 
nucleon diffusion between projectile and target, Rdiff carry information of 
’’nucleons” that have energy lower than 40MeV. Thus complimentary infor¬ 
mation is obtained from the nucleon ratios with high energy gate and the 
isospin transport ratios. 

We first perform calculations by separately varying each variable of the 
parameter set K 0 , S 0 , L , m*, fi. In the following studies, 12 parameter sets 
are constructed and listed in Table I. The average values of each variable in 
the parameter space we used are, < Ko >= 242.5 MeV, < So >= 32 MeV, 
< L >= 54.5 MeV, < m*/m >= 0.7375 and < fi >= —0.1835. Unlike 
refers. [2, 58j], there is no attempt to optimize interaction parameters to fit 
the experimental observables in this work. Deviation of < A > and < B > 
values from their ’’experimental” values is one source of uncertainties in the 
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calculated correlation coefficients Cab • As it is difficult to track the system¬ 
atic uncertainties in the transport model calculations, we do not estimate 
uncertainties in this work. 


Table 1: List of twelve parameter sets used in the IrnQMD calculations, po = 0.16/?n 3 , 


Eq = —16 MeV, and g sur = 

24.5 MeVfm 2 

i 9sur,iso — 

-4.99 MeVfm 2 

Para. 

K 0 (MeV) 

So (MeV) 

L (MeV) 

m*jrn 

fi 

1 

230 

32 

46 

0.7 

-0.238 

2 

280 

32 

46 

0.7 

-0.238 

3 

330 

32 

46 

0.7 

-0.238 

4 

230 

30 

46 

0.7 

-0.238 

5 

230 

34 

46 

0.7 

-0.238 

6 

230 

32 

60 

0.7 

-0.238 

7 

230 

32 

80 

0.7 

-0.238 

8 

230 

32 

100 

0.7 

-0.238 

9 

230 

32 

46 

0.85 

-0.238 

10 

230 

32 

46 

1.00 

-0.238 

11 

230 

32 

46 

0.7 

0.0 

12(SLy4) 

230 

32 

46 

0.7 

0.178 


To study the sensitivity of different force parameters in Skyme interac¬ 
tions on isospin observables, we calculate the covariance coefficients Cab 
based on Eq.(l)-(4) between Eve force parameters: A=K 0 , S 0 , L, rri* and 
fi and five HIC observables: B=Cl-R 2 (n/p), CI-DR(n/p), Cl-R 2 i(n/n), CI- 
R,'2 i (p/p) and Rdiff■ All the nucleon yield observables in the simulations 
are obtained for 124 Sn+ 124 Sn and 112 Sn+ 112 Sn collisions at incident energies 
of 50 and 120 MeV/u at impact parameter 2fm. For Rdiff , an additional 
mixed reaction 124 Sn+ 112 Sn is included and the calculations are performed 
with mid-peripheral impact parameter of 6 fm where a low density neck is 
formed and isospin transport between projectile and target occurs. 

Figure 1 shows the correlation coefficient Cab for Sn+Sn reactions at 
Ebeam =50 MeV/u (upper panels) and at E beam =120 MeV/u (lower panels). 
Red color bars represent positive correlations which mean observable B in¬ 
creases with parameter A, and blue color bars show negative correlations 
which mean observable B decreases with increasing parameter A. To focus 
our search for the best experimental observables which are most sensitive to 
the model parameters, we will discuss only Cab with values greater than 0.6, 
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which were represented as solid blue and red bars in Figure 1. With this 
criterion, Sq and K 0 are not very sensitive to any of the observables studied 
here. 



Figure 1: (Color online) Correlations of five obervables, Cl-R. 2 {n/p) (a), Cl-DR(n/p) (b), 
Cl-R 2 i(n/n) (c), CI-i?2i(p/p) (d), Rdiff (e) with five force parameter, Kq, Sq, L, to* and 
//. Up panels are the results for 50 MeV per nucleon, and bottom panels are for 120 MeV 
per nucleon. 


In the case of L, the slope of the symmetry energy at saturation den¬ 
sity, Cl,R diff ~ 0.7 meets this criterion. Our analysis also shows nega¬ 
tive correlations with m*, C m *,R diff ~ —0.64 to —0.70. The observation 
that Rdiff is also sensitive to m* is consistent with the results from BUU 
calculations (sol, [col [gTJ , where Rdiff values decrease with increasing m*. The 
relatively strong correlations exhibited by this observable with m* and L 
should allow one to extract the constraints of both parameters simultane¬ 
ously with reasonable uncertainties, or extract a range of L values using m* = 
0.65 - 0.9 obtained in ref. 0,0 0,0 B Ej . In ref. 0 0,0 0 Q, S 


the low energy free particle energy spectra show sensitivity to the slope of 
symmetry energy while high energy particle spectra is sensitive to the effec¬ 
tive mass splitting. When the energy cut of E c/m JA < A0MeV is applied, 
the correlation between L and CI-R^in/p): CI-DR(n/p) and CI-i? 2 i(p/p) is 
strong. This quantity will be examined in more details when the correlation 
matrix is compared to data in a future study. 

























































fi shows the strongest correlations to the experimental coalescence in¬ 
variant nucleon observables, the single and double n/p yield ratios at both 
incident energies. As shown in panels (a), (f) and (b), (g) of Figure 1, the cor¬ 
relation between the single and double n/p yield ratios and // is larger than 
0.8. The positive correlation means that the results obtained with m* < m* 
are greater than that with m* > m* which is consistent with the previous 
published results from other transport model calculationsj49|, [Ell I 52 I. [sill]. 
Strong correlation with m* is observed mainly at high energy, because the 
more violent nucleon-nucleon collisions at 120MeV/u cause larger momen¬ 
tum transfer and lead to the momentum dependent interaction (which can 
be characterized by the effective mass) to play more important roles. This 
can be understood from the expression of the density dependence of symme¬ 
try energy in Skyrrne interaction, where S(p) not only depends on isovector 
effective mass m* n but also on the isoscalar effective mass m*, i.e. S(p) = 


U F p 2/3 + A sym p + B sym p a+1 + C aym (m*, m* v )p 5/3 . 


For isoscaling ratios, the sensitivity of Cl-R 2 i(n/n) to all variables AT 0 , So, 
L, m* and fj are borderline. On the other hand, Cl-R 2 i(p/p) are particularly 
sensitive to fj due to the Coulomb repulsion as protons accelerate to higher 
kinetic energy than neutrons. Since Cl-R 2 \(p/p) depends on the difference 


of chemical potential for proton between system ”2” and ”1’ 


i.e. 


I^p2 Mpl 


which is equal to -V sym (S 2 - Si) = ~f dE k fjAS = -f dE k (^ - ^|)A 5, 
where V sym is the symmetry potential and E k is the kinetic energy of nucleon, 
it leads to positive correlation to m* and negative correlation to fj. At higher 
beam energy 120MeV/u, the influence of effective mass splitting become more 
important. It causes the strong negative correlation between Cl-R 2 i(p/p) and 
fi, and positive correlation between Cl-R 2 i(p/p) and m*. One should be able 
to extract constraints for both fj and m* from the single and double n/p 
ratios as well as the p/p ratios at 120MeV/u incident energy. Alternately, Jj 
can be extracted using the value of m*/m = 0.65 — 0.9 obtained in ref. [171. 


28|, Ll9j, 122,123, [30 


I 11 summary, by separately varying the interaction parameter sets in the 
ImQMD-Sky code, we study the influence of Ao^o^pn^and fj on the co¬ 
alescence invariant neutron and proton yield ratios at high energy region, 
for 124 Sn+ 124 Sn and 112 Sn+ 112 Sn at 50MeV/u and 120MeV/u incident en¬ 
ergy. Sensitivities to So and A'o are relatively small, Cab < 0.5. At incident 
energy of 120MeV/u, strong correlations are observed between observables 
constructed from coalescence invariant nucleon spectra with E cm jA > 40 
MeV, such as Cl -R 2 {n/p), CI-DR(n/p), and CI-A 2 i(p/p)j and the effective 
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mass splitting as well as the isoscalar effective mass, important input param¬ 
eters to the transport models. The calculations also confirm the sensitivity 
of L to the isospin diffusion observable. Since the same observable is also 
sensitive to isoscalar effective mass, this should allow one to extract the con¬ 
straints of m* and L with reasonable uncertainties. Similarly the opposite 
correlations of the nucleon yield ratios, such as CI-DR(n/p) and CI-R 2 (p/p) 
ratios to fj and m*, at 120 MeV/u reactions should allow one to disentangle 
the effects of the effective nucleon mass splitting and the isoscalar effective 
mass. 
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